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Abstract 

Understanding the volume dependence of the triton binding energy is an important step towards 
lattice simulations of light nuclei. We calculate the triton binding energy in a finite cubic box with 
periodic boundary conditions to leading order in the pionless effective field theory. Higher order 
corrections are estimated and the proper renormalization of our results is verified explicitly. We 
present results for the physical triton as well as for the pion-mass dependence of the triton spectrum 
near the "critical" pion mass, M£ nt ~ 197 MeV, where chiral effective field theory suggests that 
the nucleon-nucleon scattering lengths in the 1 So- and 3 Si-channels diverge simultaneously. An 
extension of the Liischer formula to the three-body system is implicit in our results. 
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Although Quantum Chromodynamics (QCD) is widely accepted as the underlying theory 
of strong interactions, ab initio calculations of nuclear properties in Lattice QCD remain 
one of the largest theoretical challenges in the Standard Model [H, 0] . In nuclear physics, 
the relevant degrees of freedom are pions and nucleons. Traditionally, their interactions 
are described via phenomenological potentials fitted to the nucleon-nucleon scattering data. 
More recently, the advent of model independent Effective Field Theory (EFT) approaches 
has allowed for accurate calculations of low-energy nuclear physics observables with a direct 
connection to QCD via its symmetries jsl-H]. 

In Lattice QCD, the QCD path integral is evaluated in a discretized Euclidean space-time 
using Monte Carlo simulations ■ However, this approach requires a large numerical effort 
that strongly constrains the parameters of the simulation. In particular, one is at present 
forced to use relatively small finite volumes. The energy of states calculated in the finite 
volume is shifted relative to the infinite volume limit |8j. This shift has to be taken into 
account when extracting physical observables from lattice simulations. The finite volume 
also allows for the extraction of scattering observables away from kinematic thresholds. By 
generalizing a quantum mechanical result to field theory, Liischer showed that the volume 
dependence of two-body energy states encodes the infinite volume scattering phase shift ^ 



as well as resonance properties [10J. An extension of these results to the three-body sector 
is required for the simulation of light nuclei and their scattering properties in Lattice QCD. 

This necessity motivates our investigation of the volume dependence of the triton binding 
energy within the pionless EFT. This EFT is valid for processes with typical momenta below 
the pion mass and has succesfully been applied to describe the properties of light nuclei 
[1, 0, @1- To leading order, the triton properties are determined by the nucleon-nucleon 
scattering lengths in the 1 So and 3 Si channels and a Wigner SU(4) symmetric three-body 
force . The triton has been considered in pion less EFT in a nuclear lattice formalism but 



the volume dependence was not investigated [12J. Higher-order corrections to the amplitude 



including the ones due to 2N effective range terms can be treated perturbatively I l3 | . The 
inclusion of such corrections in the pionless EFT has been studied extensively [14|-|l7|. Our 
strategy to calculate the modification of the triton in a cubic box has previously been 



employed to investigate bound states of three bosons with large scattering length [18|, [19 



In a different approach, Epelbaum and collaborators have calculated the energy of the triton 



in a finite volume by implementing a discretized version of chiral EFT on a lattice (20 
The correlation function for the three-nucleon system in the triton channel has also been 
calculated in Lattice QCD recently [2lJ , but because of the relatively large uncertainties no 
triton properties could be extracted. 

Another motivation for our work is the observation that QCD lies close to an infrared 
renormalization group limit cycle [22|. The pion mass dependence of the nucleon-nucleon 
scattering lengths from chiral EFT is compatible with simultaneously diverging nucleon- 
nucleon scattering lengths in the 1 Sq as well as in the 3 Si channel near a critical pion mass 



197 MeV [23H25| . As a consequence, it was conjectured that QCD could be tuned 



to the limit cycle by slightly changing the up and down quark masses [22[. In this scenario, 
excited states of the triton would appear near the critical quark masses. At the critical 
point itself, the triton would have infinitely many excited states with an accumulation point 
at threshold. This is a signature of the universal Efimov effect which appears in three- 



body systems with resonant short-range interactions |26|. It has recently been observed 



experimentally in atomic systems [27] and many nuclear systems can be described in an 



expansion around the ideal Efimov limit 28, |29|. The large scattering lengths make the 
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region near the critical pion mass accessible to the pionless EFT, whose input parameters 
for various pion masses can be obtained from chiral EFT. This program has success fully 



been carried out to next-to-next-to-leading order (N 2 LO) in the infinite volume 3Qj, 131 
When investigating this pion mass region with Lattice QCD, it is again desirable to have 
control over the effects of the finite volume. In particular, the influence of the finite volume 
on the infrared limit cycle is not known. 

To leading order in the pionless EFT, the Lagrangian for a system of three nucleons can 



be written as 11 



9s + 
— s A l s A 



- I W (N T T 2 a j( r 2 N) + h.c] - | [s A i (N T a 2 r A r 2 N) + h.c] + C 3 , 

where the units have been set such that h = m = 1 and <jj (t a ) are the Pauli matrices acting 
in spin (isospin) space. The degrees of freedom in this Lagrangian are the nucleon field N 
and two auxiliary dinucleon fields, tj and s A . The field tj (s A ) corresponds to two nucleons 
in the 3 Si ( 1 So)-channel. The SU(4)-invariant three-body interaction contained in £3 is 
proportional to (N^N) 3 . It can be written as a dinucleon-nucleon contact interaction with 
a dimensionless coupling constant H(A) [ill ]. 

£ 3 = - UrtforffaW + [N\t 3 a^(s A r A )N + h.c] 

V * (2) 

+ g 2 s N\s A T A y(s B r B )N 



where A is the momentum cutoff. Finite-range corrections can be incorporated as higher 



orders in the EFT [141-1171] but will not be considered here. The coupling constants g S) t 
can be matched to the two-body scattering lengths a Si t in the corresponding channel or, in 
the case of g t , to the deuteron binding energy. The three-body coupling H(A) approaches 
an ultraviolet limit cycle for large A. The analytic dependence of H on the cutoff can be 
expressed as 

cos[s log(A/A*) + arctans ] , , 

11(A) = : r, (6) 

cos [so log^A/A*) — arctan s J 

where so = 1.00624... and the phase A* has to be fixed from a three-body datum. The 
low-energy constant A* is also known as the "three-body parameter" . This renormalization 
procedure remedies the incorrect ultraviolet behavior of the EFT. Once A* has been fixed 
from a three-body datum, all other low-energy three-body observables can be predicted. For 



more details on the infinite volume case we refer the reader to the literature [111, I14J417 . 

We now consider placing the system inside a cubic box. The finite volume modifies the 
infrared regime due to momentum quantization. In a finite cubic volume with edge length L 
and periodic boundary conditions, the allowed momenta are given by ^n, n e Z 3 , yielding 
a low-momentum scale ^ corresponding to the minimal accessible non-zero momentum. 
As long as this scale is small compared to the high-energy scale given by the momentum 
cutoff A, the ultraviolet behavior of the amplitudes is unaffected by the finite volume and 
the renormalized values for the coupling constants obtained in the infinite volume can be 
used for finite volume calculations. We will explicitly verify that this claim holds for our 
results. 
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FIG. 1: Integral equation for the triton amplitudes Ft (shaded circle) and F s (shaded square). 
Single lines denote nucleons, double lines indicate dinucleons in the 3 Si channel and thick solid 
lines denote dinucleons in the 1 Sq channel. 




The properties of the triton are determined by the bound state amplitude in the spin- 1/2 
channel which is given by the integral equation in Fig. [TJ This equation has non-trivial 
solutions only for negative energies E% whose absolute values correspond to the binding 
energies. The triton amplitude has two components F t (shaded circle) and F s (shaded 
square) that correspond to the outgoing dinucleon being in the 3 Si and 1 S channels. The 
single lines denote nucleons, while the double and thick solid lines indicate full dinucleon 
propagators in the 3 Si and 1 S channels, respectively. The amplitude gets contributions from 
one-nucleon exchange as well as from the three-body interaction in £3. The loop momenta in 
this equation are quantized in a finite volume as described above. The full propagator for the 
dinucleon fields is obtained by dressing the bare propagator from Eq. ([1]) with nucleon loops 
to all orders. This yields an geometric series of diagrams that can be evaluated analytically: 



D s ,t(p ,P) 



87T 

9lt 



1 

a s ,t 



-p +p 2 /A-ie- jjj^ e 



Jez 3 
#0 



-|J1^\/-po+p 2 /4- 



(4) 



The result reproduces the infinite volume dinucleon propagator except for a volume depen- 
dent term that vanishes in the limit L — > 00. 

Using the Feynman rules from the Lagrangian in Eqs. ([[]) and fl2]) and the full dimer 
propagators from above, we can explicitly write down the coupled integral equations for the 
triton amplitude from Fig. [TJ It involves integrations over the loop energy and sums over 
the quantized loop momenta. The integrations are performed by virtue of the residue theo- 
rem while the sums over the quantized momenta are rewritten into sums of integrals using 
Poisson's resummation formula X^mez 3 ^ — m) — Xlnez 3 e 2nlx ' n , which is understood to 
be used under an integral. After projection on the quantum numbers of the triton, we have 



Hp) 



7T 



/ dv in ^ 



M 2 (y)Z{p,y) + M 3 



2H(A) 



(5) 



where Z(p, y) = [p 2 +p-y + y 2 — £3] . The matrix-valued functions M 2 (y) and M 3 (y) are 
given by 



M, 



-My) 3d « (y) 

3d t (y) -d s (y) 



M 3 (f) 



-d t (y) d s (y) 



-oL 



(6) 
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where d s ,t{y) = {gl tl^' K )F > s,t{E 3 — y 2 /2,y) depends only on the absolute value of y. The 
integral equations fl5]) obey Wigner SU(4) symmetry in the ultraviolet. As a consequence, a 
SU(4)-symmetric three-body interaction with running coupling H(A) is sufficient for renor- 
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In a finite cubic box, spherical symmetry is broken down to cubic symmetry. As a conse- 
quence, the infinitely many irreducible representations of the double cover of the rotational 
group, SU(2), become reducible in terms of the eight irreducible representations of the dou- 
ble cover of the cubic group, 2 0. The triton has j = 1/2. This partial wave is contained 



in the Gj~ representation, which also contains j = 7/2,9/2,... [32j. Assuming that the 



triton amplitude in finite volume transforms under the Gf representation, it is possible to 
decompose it into the different partial waves as 



(G+) 



(7) 



i=- 1 

J 2 ' 2 ' 



The coefficients Cum can be determined by explicitly decomposing the reducible represen- 



tations [34| and the sum over t is needed if a partial wave is contained more than once. 
Since this is not the case for partial waves less than 13/2, we will omit this index in the 
following. The vectors \jrrij) are given by \jrrij) = J2 m s C^ m i s \£U) m ) ® where the 

C's are Clebsch-Gordan coefficients, \im) is a spherical harmonic and is a spin- 1/2 
spinor. The sign in £(j) — j ± | has to be chosen such that £(j) is even in order to get the 
positive parity of the triton. 

Projecting out the Jth partial wave and performing the angular integrations in Eq. ([5]) 
yields an infinite set of coupled equations: 



F t {J \y) 



F { s J \y)j Wo 2£(J) + 



dy y 2 



M 2 (y) Z^\p, y) + M 3 (y) ^± 5 e{J) , 



F t ij) (y) 
Fs b \y), 



5jj + U'{L\n\y)< 



'( 2l(j) + l)(2f +1) 
2£{J) + 1 



E 



- 1 pi 



(n)C 



JM 



j,m+s n l(J)Q n l{J)(M-s) 



g 1 i'(M-s-m)K"'J^e^J)(M^s)^s Ly £(j) m ls C ' t(j)W '0 C 't(j)mi' '(M - s-m) 



m{t{j)),s{\) 



'JM 



(8) 

The partial waves of 



The notation m(£) is used to indicate summation over m = —£, . . . 
Z are given by 

zm(p , y) ^ Ql (l±j^i), (9) 

py V py J 

where Qt is a Legendre function of the second kind. Note that the Sjj-term in Eq. (JHJ) 
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reproduces the infinite volume result. This equation is now specialized to the case J = 1/2: 



•F^{y), 



71 



dyy 2 Yl 



M a (y) Z^(p,y) + M 3 (y)^^- 



F t u \y) 
F«\y), 



5 jh + J2^ iU) MdF\n\y) £ (- 1 ) m ^ + M^7+^ 0)m (n) 



m(<0')) 



(10) 



Since the leading term in the expansion of the Bessel functions in Eq. ([TO]) is l/(L\n\y), 
these contributions are suppressed by at least a/L. They will be small for volumes not too 
small compared to the size of the bound state. The lowest partial wave that is mixed with 
J = 1/2 is the J = 7/2 wave. Contributions from higher partial waves will be suppressed 
kinematically for shallow states with small binding momentum. This is ensured by the 
spherical harmonic in the second term of Eq. ffTUj) . Only for small lattices, i.e. when a/L is 
large, this behavior is counteracted by terms stemming from the spherical Bessel function 
je(j)(L\n\y) and higher partial waves may contribute significantly. In this first study, we will 
therefore neglect contributions from higher partial waves. We are encouraged to do so by 
our results for the bosonic case, where calculations including one more partial wave yielded 



corrections on the percent level even for small volumes [19|]. Thus, only two coupled integral 



equations remain. Details of the numerical methods used to solve these equations can also 
be found in (l9| . 

In the following, we will present our results for the energy of the triton for cubic volumes 
with various edge lengths L. In order to verify that the results are properly renormalized, 
we produced two data sets. For one set, the cutoff was set to A = 600 MeV and the three- 
body coupling was chosen such that the triton binding energy is reproduced in the infinite 
volume. For the other set, the cutoff was chosen such that the three-body force vanishes 
for the physical triton binding energy in the infinite volume. This is always possible due 
to the limit cycle behavior of the three-body force. Both data sets should agree if the 
renormalitzation has been carried out properly. 

Figure |5] shows the triton binding energy for finite cubic volumes with edge lengths L 
ranging from 17 fm down to 2 fm. The values from our two data sets (squares and circles) 
are in good agreement, indicating that our results are indeed cutoff independent. For large 
volumes, the deviation of the triton energy from its infinite volume value is small. When 
going to edge lengths smaller than about 10 fm, the energy of the state strongly decreases. 
Near L = 6 fm, the shift is already more than 100%. The dependence of the energy on L 
can be nicely fitted to a function of the form 



E 3 (L) = E 3 {L = oo) 



1 + L e -L/Lo 

L 



(11) 



indicated by the solid line in Fig. |2j We find the values Lq = 2.9 fm and c = 39 fm. The 
length scale Lq is associated with the size of the physical triton wave function. Thus the 
L-dependence of the triton energy is similar to the behavior in the two-body sector Hi. 

The triangles give the results from a lattice calculation using chiral EFT at N 2 LO by 
Epelbaum and collaborators [20]. While they agree with our results for L > 10 fm, the 
finite volume shifts are larger for smaller volumes. This discrepancy could be due to the 
fact that Epelbaum et al. calculate the volume dependence of the lowest state in the Gf 
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FIG. 2: Triton energy E3 for different edge lengths L of the cubic volume. Squares and circles are 
data sets for different cutoffs. The triangles are from a nuclear lattice calculation using chiral EFT 
at N LO [20]. The horizontal dashed line indicates the physical triton binding energy, while the 
solid line is a fit to our results (see text). 



representation, while we focus only on the j = 1/2 contribution. For small volumes, the con- 
tribution of higher partial waves is expected to be more important. The inclusion of higher 
partial waves in our approach is in principle straightforward but numerically expensive. In 
our previous work for spinless bosons [l9j], we found the contribution of higher partial waves 
to be of the order a few percent for L about three times the size of the calculated state. 
The contribution of higher partial waves in the triton will be studied in detail in a future 
publication. 

The next-to-leading order corrections to our results are given by the effective ranges in 
the 1 So and 3 Si channels. There are corrections of order r e /a and corrections of order kr e 
where r e is the effective range and k is a typical momentum. The corrections of the first 
type are of order 30%. They are dominated by the spin-triplet channel where the scattering 
length is about a factor three larger than the effective range. The corrections of the second 
type can be estimated from the typical momentum in the triton. Assuming that the binding 
momentum is shared between the nucleons in the triton, these corrections are of order 40% 
for large volumes and grow as L is decreased. Higher orders in the pionless EFT are required 
for more precise extrapolations. 

In the following, we present results for the triton and its excited states near the critical 
pion mass M% rit m 197 MeV. The pion-mass dependence of the input quantities is obtained 
from a chiral EFT calculation [30] . These input quantities are now used to determine the 
coupling constants of the "pionless" EFT, as has been described for the infinite volume case 
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31] . As already stated above, no additional input is needed to produce renormalized 
finite volume results. The pion mass dependence of the triton ground state E$ is shown in 
Fig. [3] for volumes with edge lengths L = 4.925 fm, 3.94 fm, and 2.955 fm by the upper, 
middle, and lower dashed lines, respectively. The infinite volume limit is given by the solid 
line. Again, we plot two data sets obtained using two different cutoffs as described above 
which agree well with each other. The most pronounced effect of the finite volume is a 
strong downward shift, as was expected from the results for the triton at the physical point. 
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FIG. 3: Pion-mass dependence of the triton energy E% in the vicinity of the critical pion mass 
M£ rit ~ 197 MeV. The solid line gives the infinite volume limit, while the upper, middle, and lower 
dashed line corresponds to edge lengths L = 4.925 fm, 3.94 fm, and 2.955 fm, respectively. 



Moreover, the slope of the pion-mass dependence of E$ depends on the size of the box L. 
From a linear fit to the pion-mass dependence of the triton energy in finite and infinite 
volume near the critical pion mass, we obtain the slopes AE 3 / AM n = 0.0783, 0.0467, 0.064, 
and 0.08 for the edge lengths L = oo, 4.925 fm, 3.94 fm, and 2.955 fm, respectively. As 
the volume is decreased, the slope decreases and is smallest around L = 5 fm. When L is 
decreased even further the slope increases again. 

The behavior of the first excited state when put inside a finite volume is different. For 
large volumes, the state remains again unaffected at first. When going to smaller volumes, 
the energy starts to decrease strongly. But eventually the "breakup threshold", i.e. the en- 
ergy of the relevant two-body bound state in finite volume, becomes equal to the three-body 
energy. For pion masses smaller than M° rit , the relevant two-body bound state corresponds 
to the deuteron, while for M n > M^ rit , there is a "spin-singlet deuteron" in the 1 So channel 
and the deuteron has become a virtual state. For even smaller volumes, no triton excited 
state can be found. This behavior is shown in Fig. HI The dashed lines give the energies 
of the deuteron (M n < M^ nt ) and spin-singlet deuteron (M n > M^ nt ) calculated according 
to [35||. The circles, squares, and triangles give the energies of the triton ground state, first 
excited state, and second excited state, respectively. Presumably, the triton excited state 
has crossed into the scattering regime where it would be driven away from threshold. We 
previously observed such a behavior in our investigation of three-boson bound states inside 
finite volumes [19J. Whether there is a universal relation between the binding energy in 
infinite volume and the volume size, at which the state disappears, is an interesting question 
that will be investigated in the future. The very weakly bound second excited state, that 
appears in the infinite volume for pion masses close to the critical one, could not be observed 
for the volumes with edge lengths of several fm that we investigated and, by the above rea- 
soning, has crossed the threshold already at some much larger volume. This behavior of 
the excited states makes it difficult to study the conjectured limit cycle in Lattice QCD 
calculations since very large volumes are required to observe the excited states. 

In this letter, we have studied the triton inside a finite cubic volume of edge length L with 
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FIG. 4: Pion-mass dependence of the triton spectrum near the critical pion mass for the infinite 
volume and boxes with edge lengths L = 29.6, 19.7, 14.8 fm (from left to right and top to bottom). 
The dashed lines are the energies of the deuteron (M n < M£ rit ) and spin-singlet deuteron (M n > 
M% nt ). The circles, squares, and triangles give the energies of the triton ground state, first excited 
state, and second excited state, respectively. 



periodic boundary conditions. The knowledge of the modifications of the triton energy by a 
finite volume is crucial to the understanding of lattice simulations in this channel. Using the 
framework of the pionless EFT at leading order, we have derived an infinite set of coupled 
integral equations for the partial waves of the two relevant amplitudes. This equation was 
solved for several finite volumes of an order of magnitude typical for present day Lattice 
QCD calculations. Proper renormalization of the solutions was verified explicitly. We have 
investigated the physical triton as well as the triton and its excited states near the critical 
pion mass M% nt , where the two-body scattering lengths in both channels could be tuned to 
diverge simultaneously [22] . We found that with present day volumes, the predicted excited 
states of the triton cannot be observed as bound states. 

The next step in our study is the inclusion of higher partial waves, which is straightforward 
but numerically tedious. From our previous work for spinless bosons [19], we expect the 



contribution of higher partial waves to be a few percent for L > 3a. Corrections from higher 
orders of the pionless EFT should be included for practical applications of the volume 
dependence in extrapolations. Work in these directions is in progress. So far, our calculations 
have been carried out for negative energies. While our theory is completely general, the 
numerical treatment of positive energies above the deuteron breakup threshold requires 
further work. Finally, an estimate of finite temperature effects would also be useful as 
lattice simulations are always performed at a small, non-zero temperature. 
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We note that an extension of the Liischer formula relating the infinite volume scattering 
phase shifts to the discrete energy levels in a finite volume [9| is implicitly contained in our 
work. It provides the framework to determine the low-energy constants of the pionless EFT 
in the two- and three-body sector from discrete energy levels in a cubic box. After this 
has been done, the infinite volume scattering observables can be calculated in the pionless 
EFT. An alternative is given by using a harmonic confinement instead of a cubic box. Luu 
and collaborators have shown how to reconstruct the two-body scattering phase shifts from 
the discrete energy levels in harmonic confinement 37j|. This method might be useful for 



calculations of nuclear scattering amplitudes in effective approaches with nucleon degrees of 
freedom. 

In summary, our results demonstrate that the finite volume corrections for lattice sim- 
ulations of the triton are calculable and under control. With high statistics Lattice QCD 



simulations of three-baryon systems within reach [2l[, the calculation of the structure and 



reactions of light nuclei appears feasible in the intermediate future. 
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